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Abstract 

We report a project to study scalar particles by lattice QCD sim- 
ulations. After a brief introduction of the current situation of lattice 
study of the sigma meson, we describe our numerical simulations of 
scalar mesons, a and k. We observe a low sigma mass, < m a < m p , 
for which the disconnected diagram plays an important role. For the 
ff) i kappa meson, we obtain higher mass than the experimental value, i.e. 

O \ m K ~ 2m K * 

i ! 1 Introduction 

Ph. 
D . 

r~| , The objective of our collaboration is to understand scalar mesons in the 

framework of QCD. The confidence level of the sigma meson has been in- 
creasing, and an other scalar meson k, has been reported by several 
experimental groups. In the modern hadron physics based on QCD, the 
chiral symmetry is an important ingredient and the sigma meson plays an 
essential role in it together with the pion. 

The existence of the sigma meson was obscure for many years. The 
re-analyses of the 7r-7r scattering phase shift have suggested a pole of the 
a meson with 1 = and J PC = ++ 3 . In this analysis, the chiral sym- 
metry, analyticity, unitarity and crossing symmetry are taken into account. 
Contributions of the a pole were observed in the decay processes such as 
D^ vrvrTT < and T(35) -» Ttttt 0. In 1996 PDG(Particle Data Group), 
"/o(400-1200) or a" appeared below 1 GeV mass region, and "/o(600) or a" 
in the 2002. 
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If the sigma meson exists, it is natural to consider the n meson as a mem- 
ber of the nonet scalar states of chiral SU(3)®SU(3) symmetry. Recently, 
the k with / = 1/2 is reported with mass m K ~ 800 MeV [51 [7j. 

We believe that it is very important now to investigate these scalar 
mesons by lattice QCD in order to establish scalar meson spectroscopy as 
a sound and important piece of hadron physics. Lattice QCD provides 
a first-principle approach of hadron physics, and allows us to study non- 
perturbative aspects of quark-gluon dynamics. It is a relativistic formula- 
tion, and quarks are described as Dirac fermions. It is not a model, and apart 
from numerical limitations, there are no approximations. It is not a bound 
state calculation: neither a potential model nor Bethe-Salpeter calculation. 

Lattice QCD is usually formulated in the Euclidean path integral as 

Z = J VUV^Vip e - s a-4>Dip = J vudetDe~ SG , (1) 

where Sg is the gluon kinetic action, whose continuum limit is — J d 4 xTvF^ u /4 
We construct a state with a definite quantum number and measure the decay 
in the channel, 

G{x,y) = ^ J VUV4,Vi>H(y)H( x y e~ SG -^ = (H(y)H(x)^} -> e ~ m \ x ~y\ 

(2) 

where H(x) is a hadron operator. For scalar mesons, H{x) = ip^T frfyipy 2 . 
The indices a, /j and v stand for color, flavor and Dirac indices, respectively. 

is a state whose quantum number is specified by the operator H. 
We must be very careful of the limitation of the present lattice QCD 
calculations: 

1. Sufficient statistics: Gauge configurations are generated by the Monte 
Carlo simulation, and there are statistical errors like experimental 
data. 

2. Continuum limit: Numerical simulations are performed at a finite lat- 
tice spacing, a, and we must take a — > limit at the end. 

3. Infinite volume limit: Lattice volume should be large enough to include 
a hadron. 

4. Chiral extrapolation: u and d quark masses on the current lattice are 
still large and extrapolated to zero. 

The last point may cause systematic bias, tt mesons are not sufficiently light 
and our sigma meson cannot decay into 2 tt, i.e., its width is zero. Any other 
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particle, e.g. p, A and N* , also has zero width in lattice QCD calculations 
in the literature. In the case of the sigma meson, this flaw should be kept 
in mind, since the two-pion component may be important. 

There have been several attempts at lattice study of sigma mesons. To 
our knowledge, the first such calculation was carried out by deTar and Kogut 
[S], where the so-called disconnected diagram, or the OZI forbidden type 
diagram is discarded. The channel was called 'valence' sigma, ay. They 
measured the screening masses and observed that ay is much heavier than 
the 7r meson at the zero temperature, while ay and ir degenerate as the 
temperature increases over T c . Kim and Ohta calculated the valence sigma 
mass with staggered fermions for the lattice spacing a = 0.054 fm and lattice 
size 48a = 2.6 fm [S]. They obtained m^/m^ = 1.4 ~ 1.6 by varying 
m n /m p = 0.65 ~ 0.3. 

Lee and Weingarten ^U] have stressed the importance of the mixing the 
scalar meson and glueballs and concluded that /o(1710) is the lightest scalar 
glueball dominant particle, while /o(1390) is composed of mainly the u and 
d quarkonium. Alford and Jaffe analyzed the possibility that the sigma 
particle is an exotic state, i.e., qqqq by a lattice QCD calculation [TT] . 

All these calculations are in the quench approximation, i.e., the fermion 
determinant in Eq.Q is dropped, which corresponds to ignoring quark pair 
creation and annihilation diagrams. McNeile and Michael observed that the 
a meson masses in the quench approximation and in the full QCD simulation 
are very different |12j . They also considered the mixing with the glueballs. 
They obtained a very small sigma meson mass, even smaller than the ir 
mass, in the full QCD case. 

There are two ongoing projects of a meson spectroscopy : Riken-Brook 
haven- Columbia (RBC) collaboration jOj and Scalar collaboration |14l I15j . 
The two approaches are complementary. The RBC collaboration employs 
the domain wall fermions, which respect the chiral symmetry, but include 
a quench approximation, while Scalar collaboration uses Wilson fermions, 
which break the chiral symmetry at a finite lattice spacing, but performs 
the full QCD simulation. RBC reported their simulation at a -1 = 1.3 GeV 
on a 16 3 x 32 lattice. They remedied the quench defect with the help of the 
chiral perturbation. They observed that the masses of the non-isosinglet 
scalar (ao) and the singlet scalar (a) are almost degenerate when the quark 
mass is heavy (above s quark mass regions), and that as the quark mass 
decreases, the ao mass remains almost constant, but the mass of the a 
decreases. 



3 



2 Propagators 



The quantum numbers of the a meson are 1 = and J PC = ++ ; we adopt 
the a meson operator of 

„. , u(x)u(x) + d(x)d(x) . 
a(x) = -j= , (3) 

where u and d indicate the corresponding quark spinors, and we suppress 
the color and Dirac indices. The a meson propagator is written as 

G a (y,x) = | J ' VUVuVuVdVd[a{y)a{x)^)e- SG - nDu - dDd . (4) 

By integrating over u, u, d and d fields, the a meson propagator is given by 

G a (y,x) = - (TrD- 1 (x,y)D- 1 (y,x)) 

+2((a(y) - (a(y)))(a(x) - (a(x)))} (5) 

where <r(x) = TrD~ 1 (x, x). "Tr" represents a summation over color and 
Dirac spinor indices. In Eq.Q, -D -1 's are u and d quark propagators. Here 
we assumed that the u and d quark propagators are equivalent because u 
and d quark masses are almost the same. From Eq.(|5jl. we can see that 
the a propagator consists of two terms. The first term corresponds to the 
connected diagram, i.e., a qq-type meson. 

The second term is the "disconnected" diagram; it is the correlation of 
a = TvD^ 1 at two points x and y. The term "disconnected" is not appro- 
priate since the corresponding matrix element is, of course, not factorized; 
quark lines are connected by gluon interactions. Nevertheless, we use this 
jargon in the following. 

The quantum number of the a meson (1 = and J p = 0+) is the same 
as that of the vacuum, and the vacuum expectation value of the a operator, 
(a(x)), does not vanish. Therefore, the contribution of (<t(x)) should be 
subtracted from the a operator. 

For k + , the operator is H(x) = J2 a J2u^u u u^ an d we have only the 
connected diagram, 

G K (y,x) = -{TYD-\x,y)D-\y,x)), (6) 

where D~ l is the s quark propagator. 



4 



3 Numerical simulations 



In this project, we intend to use standard well-established techniques for 
numerical calculations, and want to see the outcome. We employ Wilson 
fermions and the plaquette gauge action. Full QCD simulations were done 
by the Hybrid Monte Carlo (HMC) algorithm. 

CP-PACS performed a very large-scale simulation of light meson spec- 
troscopy in the full QCD calculation ^B]. We use here the same values of the 
simulation parameters, i.e., (3 = 4.8 and k = 0.1846, 0.1874, 0.1891, except 
lattice size; our lattice, 8 3 x 16, is smaller. We employ the point source and 
sink; smaller lattice size leads to larger mass due to a higher state mixture. 
In other words, our mass values on the small size lattice should be consid- 
ered as the upper limit. We have checked that the values of m n and m p are 
consistent with those of CP-PACS. 




Figure 1: Propagators of tt, p and a for n = 0.1891. 

It is very difficult to evaluate the disconnected part of the propagator, 
since we must calculate TrD~ 1 (x, x) for all lattice sites x. We used the 
noise method to calculate the disconnected diagrams and the subtraction 
terms of the vacuum (a). Each of these terms is of the order of ten, and 
{(a — (<t))(<t — {&))} becomes less than 10 -4 , as shown in FigEJ Therefore, 
high accuracy is required for the calculation. One thousand random Z2 
numbers are generated. Our numerical results show that the values of the 
first and the second terms in Eq.@ are of the same order. Therefore, in 
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order to obtain the signal correctly as the difference between these terms, 
high-precision numerical simulations and careful analyses are required. We 
have investigated the relationship between the amount of Z2 noise and the 
achieved accuracy in Ref. |14| . Gauge configurations were created by HMC 
in the SX5 vector supercomputer, and most disconnected propagator cal- 
culations by the Z2 noise method were mainly performed on the SR8000 
parallel machine at KEK. 

The propagators of tt, p and a for ft = 0.1891 are shown in Fig^ The 
connected and disconnected parts of the a propagator are shown in FigEl 
It is difficult to obtain a propagator at large r since the precision of our 
calculation is limited to 0(G(t)) ~ 10~ 4 . 
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Figure 2: Propagators of the connected and disconnected diagram of a for 
k = 0.1891. 

From our results, we evaluate the critical value of the hopping parameter 
k c = 0.195(3) and lattice space a = 0.207(9) fm (CP-PACS has obtained 
k c = 0.19286(14) and a = 0.197(2) fm). Figure ED (left) shows masses of 
p, (j and 7T as functions of 1/k. We find that m a /m p at the chiral limit is 
0.33 ± 0.09. 

Finally, we present our preliminary result for k meson using the same 
common configurations. The s quark is treated as a valence, i.e., is used 
only in the propagator Q, not quark. We adopt the same hopping 

parameter values, k = 0.1846, 0.1874 and 0.1891 for u and d quarks. We 
calculate three values of the hopping parameter for the s quark: k s = 0.1835, 
0.1840 and 0.1845. For each k s , we calculate masses of k, K* and K mesons, 
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and extrapolate them to the chiral limit. Then we tune the s quark hopping 
parameter, k s , to give the best experimental values for mx* and mx- Figure 
01 (right) shows m K a, rriK'O- and m^a as functions of 1/k for k s = 0.1840. 
Our preliminary analysis shows that value of m K /mK* at the chiral limit is 
around 2.0. 
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Figure 3: Left : m„, m a and in the lattice unit as functions of the 
inverse hopping parameter. Right : m K , tjik* and mx in the lattice unit as 
functions of the inverse hopping parameter. The s quark hopping parameter 
is k s = 0.1840. 



4 Concluding Remarks 

We have reported our exploratory study of the scalar mesons based on the 
full QCD lattice simulation with dynamical fermions. Our results indicate 
the existence of a light a in the region m n < m a < m p . 

An interesting observation is that the disconnected part gives a signifi- 
cant contribution, and this diagram makes the a meson light. This cannot 
be accessed in the framework of the non-relativistic quark model. This point 
should be kept in mind in future phenomenological analyses of the sigma 
meson, and also in lattice studies. Note that the k meson does not have such 
a mechanism, and therefore these two scalar meson masses can be different. 
The n and the valence a (connected part) have the same structure in their 
propagators, but s quark mass is heavier than those of u and d. Therefore 
m K > m conn , where m conn is a mass corresponding to the connected part, 

The calculations reported here have limitations discussed in Sect.l. We 
expect that these defects will be gradually overcome. In lattice calculation, 
once someone establishes the scale of simulation required to obtain mean- 
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ingful results, many improved and large-scale works succeed. Therefore in 
a few years, the lattice study of scalar mesons will provide important and 
fundamental information to deepen our understanding of hadron physics. 
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